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1.  INTRODUCTION 


In  this  report,  we  will  derive  the  equations  of  motion  for  a 
single-degree-of-freedom,  rate-integrating  gyro,  taking  into  account 
asymmetry  of  design,  manufacturing  tolerances,  and  assembly  tolerances 
which  result  in  non-ideal  geometries  and  dynamics.  We  will  also  review 
the  test  methods  and  the  data  reduction  procedures  used  in  determining  the 
performance  coefficients  of  a  gyro.  Lastly,  we  will  briefly  investigate  the 
torques  due  to  convection  currents. 

There  are  two  complementary  procedures  commonly  used  for 
determining  the  performance  coefficients  which  we  will  review;  they  are: 

a)  tumble  tests  (either  continuous  rotation  or  discrete  positions)  and 

b)  precision  linear  vibration  tests.  In  theory,  a  precision  centrifuge  with 
a  counter-rotating  table  (not  necessarily  a  1:1  counter-rotation)  may  be 
used.  However,  in  practice,  the  results  are  extremely  sensitive  to  align¬ 
ment  of  the  counter-rotating  table  axis  to  the  centrifuge  axis,  and  therefore, 
that  procedure  will  not  be  pursued  further. 

Data  from  tumble  tests  and  linear  vibration  tests  are  commonly 
reduced  by  using  formulae  which  are  partially  or  fully  computer-automated. 
There  is  seldom  any  thought  given  to  the  formulae  used  nor  their  applica¬ 
bility  to  a  particular  gyro  design  or  its  application.  The  simplifying 
assumptions  (there  are  always  such  assumptions)  used  in  the  data  reduction 
procedures  may  or  may  not  be  appropriate  for  the  particular  design  or  the 
application. 

One  of  the  principal  purposes  of  this  report  is  to  review  the 
assumptions,  explicit  or  implicit,  in  the  choice  of  test  method  and/or  data 
reduction  procedure.  The  test  'ngineer  must  then  review  the  error  sources 
and  use  his  judgment  as  to  the  appropriateness  of  the  test  methods  and  the 
data  reduction.  In  some  cases,  analyses  and/or  tests  may  be  required 
prior  to  making  a  final  judgment. 
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2.  MISALIGNMENTS  WITHIN  A  GYRO 


2.  1  Introduction 

In  this  section,  we  will  investigate  the  effects  of  non-ideal 
geometries  arising  from  design  compromises,  material  variations,  manu¬ 
facturing  tolerances,  and  assembly  tolerances.  It  is  usually  assumed  that 
the  true  spin  axis  (SA)  as  determined  by  the  bearing  geometry,  is  normal 
to  the  true  output  axis  (OA)  as  determined  by  the  gimbal  pivot  or  suspension 
geometry.  The  input  axis  (IA)  is  then  taken  to  be  normal  to  SA  and  OA. 

Thus,  the  three  axes  form  an  orthogonal  coordinate  system.  It  is  also 
generally  assumed  that  IA,  OA,  and  SA  are  principal  axes  of  inertia  for 
the  float.  In  reality,  none  of  these  assumptions  is  likely  to  be  exactly  true. 

2.  2  Coordinate  Systems 

There  are  several  coordinate  systems  that  could  be  used  in  this 
study,  each  having  its  relative  advantages  and  disadvantages.  We  will  use 
three  basic  coordinate  systems  in  our  development:  a)  the  reference  sys¬ 
tem  attached  to  the  gyro  case  as  indicated  by  reference  markings  and/or 
reference  surfaces,  b)  one  attached  to  the  gyro  float,  and  c)  an  inertially- 
fixed  coordinate  system.  Other  auxiliary  coordinate  systems  will  be  used, 
as  convenient,  such  as  an  earth-fixed  coordinate  system. 

In  any  single-degree-of-freedom  gyro,  there  are  two  fairly 
well  defined  axes;  they  are  the  spin  axis  SA  and  the  output  (pivot)  axis  OA. 
Though  these  axes  may  be  well  defined  by  geometry,  they  are  not  necessarily 
orthogonal,  nor  are  they  easy  to  locate  in  an  assembled  gyro.  Assuming 
that  the  pivot  clearances  are  very  small,  as  is  usual,  then  any  angular 
motion  of  the  gimbal  relative  to  the  gyro  case  must  be  about  OA  only, 
whether  or  not  OA  is  perpendicular  to  SA. 

The  three  reference  axes,  which  are  indicated  by  reference 
markings  and/or  reference  (mounting)  surfaces  on  the  gyro  case,  are: 
a)  the  input  reference  axis  (IRA),  b)  the  output  reference  axis  (ORA),  and 
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c)  the  spin  reference  axis  (SRA).  These  are  nominally  parallel  to  LA,  OA, 
and  SA,  respectively,  and  form  an  orthogonal  coordinate  system  such  that 
the  vector  cross  product 

T  x  J  =  E  (2-1) 

where  T,  j”,  and  Tc  are  unit  vectors  along  IRA,  ORA,  and  SRA,  respectively. 
The  origin  of  this  coordinate  system  is  taken  to  be  at  the  centroid  (C)  of  the 
gimbal  (float)  when  at  electrical  null  and  at  neutral  buoyancy.  For  con¬ 
venience,  IRA,  ORA,  and  SRA  will  be  designated  1,  O,  and  S,  respectively, 
as  shown  in  Fig.  2-1. 

For  the  gimbal,  we  will  choose  an  XYZ  coordinate  system  with 
the  origin  at  C,  the  Y  axis  parallel  to  the  true  pivot  axis  (OA)  and  the  YZ 
plane  parallel  to  the  true  spin  axis  (SA).  Let  C  rx  be  the  angle  between 
SA'and  the  Z  axis,  as  indicated  in  Figure  2-1,  where  SA'  is  the  projection 
of  SA  on  the  YZ  plane. 

We  will  define  the  rotor  as  being  the  gyro  wheel  plus  any  asso¬ 
ciated  rotating  parts  such  as  bearing  races  and  ball  retainers.  The  gimbal 
or  float  is  defined  as  the  complete  assembly,  including  the  rotor,  that 
rotates  as  a  unit  about  OA  relative  to  the  gyro  case. 


i 


Fig.  2-1.  Gyro  Reference  Axes  and  Gimbal  Coordinate  System 
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Moments  of  Inertia 


2.  3 

Moments  of  inertia  and  products  of  inertia  arise  when  we  apply 
Newton's  Second  Law  of  Motion  to  a  rotating  rigid  body.  It  is  assumed 
that  the  reader  is  familiar  with  the  definitions  of  moments  of  inertia  and 
products  of  inertia,  and  we  will  review  only  those  additional  facts  which 
are  pertinent  to  our  study  of  the  gyro. 

For  any  rigid  body,  with  or  without  axes  of  symmetry,  there 
are  always  at  least  three  orthogonal  axes  passing  through  the  centroid  for 
which  the  products  of  inertia  vanish.  These  axes  are  called  principal  axes 
and  the  moments  of  inertia  about  these  axes  are  known  as  principal  moments 
of  inertia.  Any  axes  of  symmetry  would  be  principal  axes  since  the  pro¬ 
ducts  of  inertia  about  those  axes  would  vanish. 


The  moment  of  inertia  of  the  gimbal  about  any  line  p  passing 

through  the  centroid  may  be  given  in  terms  of  the  principal  moments  of 

inertia  J  ,  J  ,  and  J  about  three  orthogonal  principal  axes  U,  V 

and  W  and  the  direction  cosines,  1  ,  i  ,  and  £  of  the  line  p  with 

pu  pv  pw 

respect  to  the  principal  axes. 


J  =  J  i2  +  J  Jt2  +  J  i2 
p  uu  pu  w  pv  WW  pw 


(2-2) 


The  principal  moments  of  inertia  are  stationary  values  of  Eq.  (2-2),  i.  e. , 
Jp  is  a  maximum  or  a  minimum  when  the  line  p  coincides  with  a  principal 
axis. 


Assume  that  wc  have  determined,  by  some  means,  the  principal 
moments  of  inertia  of  the  gimbal,  and  that  the  principal  axes  U,  V,  and  W 
arc  nearly  parallel  to  the  1,  C,  and  S  axes,  respectively ,  a:  is  usually  the 
case.  Let  9  ,  t  ,  and  6  be  the  successive  small  Euler  antic  rotations 

u  v  V 

about  the  U,  V,  and  W  axes,  respectively,  required  to  rotate  the  UVW 
coordinate  system  into  the  lOS  coordinate  system,  it  can  be  shown  that 
the  direction  cosines  are? 


(2-3) 


*Iu  ~  1  Ov  “  *Sw 


l*.  ~-tc  *-  e 

Ow  Sv  u 


-  -iT  =-0 
Su  Iw  V 


*jv  *■  -  W  *  e. 


Ou  °w 


The  three  principal  moments  of  inertia  of  the  gimbal  usually  dif- 

fer  from  one  another  by  less  than  an  order  of  magnitude  and  the  principal 

axes,  in  general,  are  nearly  coincident  with  the  I,  O,  and  S  axes,  If6u» 

6  ,  and  G  are  each  less  than  0, 05  rad  (<  3°),  the  moments  of  inertia  of 
v  w 

the  gimbal  about  the  reference  axes,  as  determined  from  Eqs.  (2-2)  and 
(2-3),  are  approximated,  without  significant  error,  by: 

j.-  j  +  j  e2  +  j  e2  -  j 

1  uu  W  W  WW  V  uu 


j  ~ j  e2+j  +j  e2 -  j 

O  uu  W  W  WW  °u  w 


(2-4) 


jc  -  j  e2  +  j  e2  +  j 

S  uu  V  W  U  WW 


WW 


The  above  results  should  not  be  surprising,  since  the  principal 
moments  of  inertia  are  stationary  values  with  zero  slope. 
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3.  DYNAMICS  OF  THE  GIMBAL 


3.  1  Introduction 


Any  system  of  force*  and  torques  applied  to  a  rigid  body  may  be 

replaced  by  a  single  resultant  force  7  acting  through  the  centroid  and  a 

c 

single  resultant  torque  Tc«  For  example,  consider  the  single  force  7^ 
shown  in  Fig.  3-  1(a)  which  is  eccentric  to  the  centroid  C. 


Fig.  3-1.  Eccentric  Applied  Force 

Apply  two  opposed  forces  and  7^  at  C  in  a  line  parallel  to 

the  line  of  action  of  ^  as  shown  in  Fig.  3-l(b)  and  such  that  |^jc|  = 

|7jcJ  =  |Fj|.  Since  the  vector  sum  of  Tjc  and  Fjc  is  zero,  the  net  force 

acting  on  the  body  is  unchanged.  However,  we  can  now  consider  the  force 

system  to  be  a  force  acting  through  the  centroid  and  a  torque  7jh  due 

to  the  couple  formed  by  the  for  :es  and  The  same  procedure  may 

be  followed  for  each  non-central  force.  All  the  resultant  forces  acting 

through  the  centroid  may  be  summed  vectorially  to  obtain  but  a  single 

resultant  force  ^  and  all  the  torques  may  be  summed  vectorially  to  obtain 
c 

a  single  resultant  torque  Tg. 

The  resulting  motion  of  the  rigid  body  may  be  considered  as  a 
linear  translation  of  the  centroid  due  to  the  resultant  force  and  a  change  in 
angular  momentum  due  to  the  resultant  torque  T  .  Assuming  a  constant 
mass,  the  differential  equations  of  motion  are: 


7 


Ci-l) 


F  =  d(mV  )/dt  =  ma 
c  c  c 

Tc  a  dJT/dt  (3-2) 

whore 

m  a  mui 

V  *  velocity  vector  o £  centroid 

T  a  acceleration  vector  of  centroid 
c 

IT  a  angular  momentum  vector 
t  *  time 

Equation  (3-1)  atates  that  the  rate  of  change  of  the  linear  mo¬ 
mentum  (mVc)  is  equal  in  magnitude  and  direction  to  the  resultant  applied 

force  Y  acting  through  the  centroid.  Equation  (3-2)  states  that  the  rate 
c 

of  change  of  the  angular  momentum  (IT)  is  equal  in  magnitude  and  direction 
to  the  resultant  applied  torque  T^.  Equation  (3-2)  is  easily  derived  from 
the  basic  Eq.  (3-1). 

Though  we  are  primarily  interested  in  he  angular  motions  of 
the  gimbal  and  will  stress  the  use  of  Eq.  (3-2).  tl».  gimbal  structure  is  not 
a  perfectly  rigid  body  and  compliance  torques  resulting  from  linear  acceler¬ 
ation  of  the  centroid  will  also  be  considered.  Of  course,  a  gyro  designer 
must  consider  the  translational  forces  acting  on  the  pivots,  rotor  bearings, 
and  other  parte  for  strength  as  well  as  for  minimizing  deflections. 

3. 2  Principal  Axes  of  Rotor  and  Gimbal  Structure 

Because  of  tight  manufacturing  tolerances,  accurate  dynamic 
balancing,  and  extreme  care  in  assembly,  the  spin  axis  (SA)  is  a  principal 
axis  of  the  rotor.  The  other  principal  axes  of  the  rotor  must  be  in  a  plane 
normal  to  SA  and  passing  through  the  centroid  of  the  rotor.  We  will  as¬ 
sume  that  all  axes  in  that  plane  and  passing  through  the  centroid  are 
principal  axes  of  inertia  and,  therefore,  the  moments  of  inertia  about 
these  axes  are  all  equal.  The  above  assumption  is  true  for  the  vast 


•  8* 


majority  of  gyros.  There  are  a  few  gyros,  such  as  the  "Nutatron",  which 
are  purposely  designed  to  have  a  rotor  with  unequal  moments  of  inertia 
about  their  two  radial  principal  axes  of  inertia,  but  such  gyros  are  not 
yet  commercially  significant. 

Manufacturing  and  assembly  tolerances  for  the  gimbal  do  not 
permit  as  accurate  a  knowledge  of  the  location  of  the  principal  axes  as  we 
had  for  the  rotor  alone.  However,  gimbals  are  generally  designed  so  that 
the  principal  axes  are  nominally  parallel  to  the  X,  Y,  and  Z  axes  and 
the  Y  axis  passes  through  or  close  to  the  centroid  C  of  the  gimbal.  There¬ 
fore,  we  may  use  Eqs.  (2-4)  for  the  gimbal,  without  appreciable  error,  for 
the  vast  majority  of  gyros  and  we  will  so  assume. 


3.  3  Angular  Momentum  of  Rotor  and  Gimbal 

Since  the  only  axes  externally  available  are  the  reference  axes 
(I,  O,  and  S  in  Fig,  2-1),  we  will  determine  the  angular  momentum  of 
the  gyro  with  respect  to  these  axes.  Consider  first  the  angular  momentum 
of  the  rotor  relative  to  the  gimbal. 

The  angular  momentum  components  of  the  rotor  relative  to  the 
gimbal  in  the  XYZ  coordinate  system  are: 


where 


H  =  0 

rx 


H 


rx 


•J  Uj  * 

r  r  rx 


-H  ; 
r  rx 


H  *  J  u.  =  H 
rz  r  r  r 


(3-3) 


H  =  J  it  =  angular  momentum  of  the  rotor  (plus  asso- 
r  r  dated  rotating  parts)  relative  to  the  gimbal 

J  *  moment  of  inertia  of  the  rotor  about  its  axi6  of 

x* 

rotation 

uj r  =  angular  velocity  of  the  rotor  relative  to  the  gimbal 
£  =  misalignment  angle  as  defined  in  Fig.  2-) 


In  order  to  determine  the  angular  momentum  components  of  the 
rotor  relative  to  the  gimbal  in  the  IOS  system,  we  will  perform  the  follow¬ 
ing  transformations: 

a)  Rotate  the  XYZ  system  about  the  Y  axis  through  the  small 
anglo  -cpq  required  to  bring  the  gimbal  back  to  its  electrical 
null  position, 

b)  Rotate  the  resulting  system  successively  through  the  small 
Euler  angles  C  ,  Cy,  and  £.  to  bring  it  into  coincidence  with 
the  ICS  system  which  is  fix&d  to  the  gyro  case. 

Since  the  rotation  angles  are  all  small  («  1  rad),  we  may  use 
small  angle  approximations  for  the  sines  and  cosines  and  neglect  products 
of  the  sines.  We  find  the  angular  momentum  of  the  rotor  in  the  IOS 
system  relative  to  the  gimbal  is: 

Rr  “  H,<»o*  Cy|T  *  Hr(Cx  ‘  ‘rx1  T  +  Hr^  <3-‘U 

Assuming  the  principal  axes  of  inertia  of  the  gimbal  are  nominal¬ 
ly  parallel  to  the  I,  O,  and  S  axes,  then  using  Eqs.  (2-4),  the  angular 
momentum  of  the  gimbal  with  the  rotor  locked  in  the  IOS  system  is: 

^g  Ji  tt‘i 1  f  Jo^tt*o  +  ^O^  J  +  JS  “S*  0-5) 

where 


■j»j,  u;q,  uog  =  angular  velocity  components  of  the  gyro  case 
relative  to  inertial  space  about  the  I,  O,  and  S 
axes,  respectively  (see  Fig.  2-1). 

<p_  =  angular  velocity  of  the  gimbal  relative  to  the  gyro  case 
°  about  OA  and  «■  1. 

The  total  angular  momentum  of  the  gimbal  is  the  sum  of  Eqs. 
(3-4)  and  (3-5).  As  a  mnemonic  device,  with  no  loss  in  generality,  let 
(^x  •  Crx)  -  Cj  &od  C  =  C  q»  then  the  total  angular  momentum  of  the  gimbal 
is: 


—  1 0< 


(3-6) 


TT  «  TTr  f  TTg  -  [Hr(<p0  -  cQ)  +  Jx  «j]T 

+  [h.  Cj  +  J0(«>0  +  Vo' X  +  ^Hr  +  Js  ^ 


3. 4  Torque  Equations 

The  torque  components  about  the  I,  O,  and  S  axes  required  to 
obtain  changes  in  angular  momentum  are  determined  by  applying  Eqs.  (3-2) 
and  (3-7),  below,  to  Eq.  (3-6). 

di/dt  =  tt)sf  -  <BqTT 

dj/dt  =  u>jlc  -  <^r  (3-7) 

dlc/dt  =  U.J  -  ulJ 

TI  ~  Hr<“>0  +  *0  ~  ‘i  %)  +  Jx  •«>!  +  Jr  ir(^0  -  C0> 

"  Jo*o  “s  •  (Jo  •  Js)  “o  “s  (3-8) 

T0  ^  “Hr  ^“i  *  +  J0^“,0  +  *^0^ 

+  Jr  Cx  ir  -  (Js  -  Jj)  u»s  ««j  (3-9) 

r  i  *  * 

Tg  -  -Hf  [(Vq  -  Cq)  tt>Q  -Cj«^J^JgWgfJr 

+  Jo*o  "i  *  (JI  •  Jo)  “i  *0  {i’l0) 

It  i3  customary,  for  purposes  of  comparison,  to  normalize  the 
torque  components  by  dividing  them  by  the  angular  momentum  H  of  the 
rotor  relative  to  the  gimbal.  It  is  also  common  practice  to  break  out  the 
normalized  viscous  drag  torque  (-C^q/H^,)  which  opposes  the  angular 
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velocity  cpQ,  the  normalized  elastic  restraint  torque  (-Ke  Pq/H  )  which 
oppose 8  the  angular  displacement  cp^,  and  lastly,  the  normalized  command 
torque  (-Kji)  where: 

C  =  damping  coefficient 

K  =  elastic  restraint  coefficient 
e 

Kt=  torquer  scale  factor 

i  =  torquing  current 

/ 

4 

The  sign  convention  adopted  for  the  torquing  current  is  such 
that  a  positive  current  will  exert  a  negative  torque  about  OA  or,  equivalently, 
will  cause  a  posit  we  precession  rate/about  IA  when  the  gimbal  is  maintained 
at  its  null  position.  The  normalized  torques  about  the  1,  O,  and  S  axes  are: 

• 

TT  *  M  =(wQ  t  <PC  -  Cj  U)g)  +YP  <Sj  *  (i'o  -  Cq) 
r  r  r 

T  i  -  T 

Ol  ‘O  JS  ,,  ,14 

-w* o  *s  *  —  “o  (3‘U) 


=  Mo  ’  IT  *0  "  IT^o  -kt1  =  "^i  "  ((po  "  co)  *s  ] 

•  ••  tt,r  JS  “  Jj 

+tt,(“o  +cpo)  +  ci « - 1 “fc  "i  (3'12) 

r  r  r 


TT  =  M s  =  4<<po  -  co>  “o  ■  CI  “iJ  +  if  “s  *  T- 
r  r  r 


O  * 

frv  0*1” 

r  i 


JI  "  J0 

R  "to 


(3-13) 


where  Mj  and  Mg  are  the  total  applied  normalized  torque  components  about 
the  I  and  S  axes,  respectively,  and  MQ  is  the  total  normalized  torque  about 
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the  O  axis  plus  the  normalized  viscous  drag  torque*  the  normalized  elastic 
restraint  torque*  and  the  normalized  command  torque*  i.  e,  * 


(3-14) 


If  the  misalignment  angles  Cj  end  Cq  are  negligible  and  if  the 
gimbal  is  maintained  at  its  null  position  (<Pq  =  0)*  then  Eqs.  (3-11)*  (3-1?.)* 
(3-13)  become  the  following  well-known  equations: 


w  _  .  Ji  *  Jo-Js 

MI  =  mO  +  TT  “I  '  “TT—  “O  “fe 


M0-KTi  »  -  VT7  iQ  -  “V"1  “fc  "I 

e 

Jc  •  "r  JT  “ 

Ms  =  ir  -s *  it  -  -hr-2  “i  -o 

r  v  r 


(3-15) 


Note  that  if  «»>q  -  s»g  =  i  =  0,  then  from  the  second  of  Eqs.  (3-15),  we  see 
that  a  positive  applied  normalized  moment  Mq  will  result  in  a  negative  pre¬ 
cession  rate  when  the  gimbal  is  maintained  at  its  null  position. 
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4.  MODEL  EQUATION 


In  most  application*  and  testa,  the  gimbal  is  maintained  at  or 

close  to  its  electrical  null  position  by  a  platform  or  table  slaved  to  the  gyro 

and/or  by  a  feedback  torquing  current.  The  null  error  will  depend  on  the 

* 

tightness  of  the  servo  loop.  We  may  rewrite  Eq.  (3*12)  with  *  cp^  =  0 
and  Mq  expressed  in  black  box  form  as  a  function  of  acceleration.  Varia¬ 
bles  other  than  acceleration,  such  as  power  supply  variations  and  tempera¬ 
ture  variations,  are  not  usually  included  in  the  model  equation,  but  are 
treated  as  random  variables  which  result  in  uncertainties  in  the  drift  rates. 
Unfortunately,  there  are  several  model  equations  in  common  use.  Some  are 
expressed  in  terms  of  drift  rates,  others  are  expressed  in  terms  of  normal¬ 
ized  torques  which  have  the  dimensions  of  drift  rates  but  are  of  opposite 
sign.  Some  equations  are  expressed  in  terms  of  acceleration  inputs,  where¬ 
as  others  are  expressed  in  terms  of  specific  force  inputs  which  have  opposite 
signs.  Therefore,  it  is  hardly  surprising  that  there  is  some  confusion,  and 
that  some  model  equations  are  incorrect. 

Up  to  this  point,  we  have  not  specified  any  dimensional  units  in 

the  various  equations,  but  merely  expected  them  to  be  consistent.  For 

example,  Eqs.  (3-11),  (3-12),  and  (3-13)  would  be  dimensionally  consistent 

if  the  applied  torques  are  in  dyne-cm,  the  angular  momentum  in  gm  •  cm4/ 

2 

a,  the  moments  of  inertia  in  gm  •  cm  ,  the  angles  in  radians,  the  angular 

2 

rates  in  rad/s,  and  the  angular  accelerations  in  rad/s  .  In  our  model 
equation,  we  will  adopt  the  more  common  practice  of  expressing  u>j,  a  q, 
and  u)  g  in  deg/hr  and  <«q  in  (deg/hr)/s.  We  will  leave  u>r  and  u)f  in  rad/ 
s  and  rad/s^.  We  will  also  adopt  the  practice  of  expressing  the  acceler¬ 
ation  in  units  of  the  local  g  value.  The  model  equation  and  the  definition  of 
symbols  follow. 

®I  =  °F  +  DI*I  +  DOaO  +  DS*S  +  DII*I  +  DOOaO  +  DSSaS 


+  DIO  *1*0  +  D0S  aO*S  +  DIS  *I*S 
Jq  •  Wr  ^ S  ~ 

+  !T“  “o  +  206'  300  CI  - H — 

r  r  r 


+  KTi  -  C 

1067306 


o 


tt). 


(4-1) 


Precriing  page  link 
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Note  that  the  1  deg/hr  =  1  aec/a  and  1  rad/a  206,300  aec/a* 

in.,  u>q,  a»s  =  angular  rate  componenta  of  the  gyro  caae  about 
A  the  reference  axes  with  reapect  to  inertial 

apace  -  deg/hr 

d>n  =  angular  acceleration  of  the  gyro  caae  about 

u  ORA  -(deg/hr)/a 

aj,  aQ,  ag  =  componenta  of  total  applied  acceleration,  a, 

along  the  reference  axea  -  £ 

a  s  aT  -  £  =  applied  acceleration  vector  -  £ 

a_  =  total  acceleration  vector  with  reapect  to  inertial 

apace  -  £ 

£  -  gravity  vector  (in  direction  of  maximum  gravity 

gradient)  -  £ 

£  =  unit  of  acceleration  equal  in  magnitude  to  the 

local  (or  atandard)  gravity  value 

Dp  =  fixed  reatraint  (bias)  drift  rate  -  deg/hr 

D.,  Dn,  D-  =  linear,  acceleration-sensitive  drift  rate  coef- 

ficienta  -  (deg/hr)/£ 

DrT,  D__,  Dcc  =  quadratic,  acceleration-sensitive  drift  rate 

n  coefficient  -  (deg/hr)/£2 

D.  Doc,  Dtc  =  cross -coupled,  acceleration-sensitive  drift  rate 

*  coefficients  -  (deg/hr)/£^ 

K,j.  =  command  rate  scale  factor  (deg/hr)/A 

i  =  command  rate  current  -  A 

J.,  Jq,  Jg  =  momenta  of  inertia  of  the  gimbal  about  the  refer¬ 

ence  axea  -  gm  •  cm2 

a)  =  angular  rate  of  the  rotor  relative  to  the  gimbal  - 

r  rad/a 

e 

u>  =  angular  acceleration  of  the  rotor  relative  to  the 

gimbal  -  rad/a2 

H  =  angular  momentum  of  the  rotor  relative  to  the 

r  gimbai  -  gm  •  cm*/. 
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ci»  co 


misalignment  of  SA  with  respect  to  SRA  about  IRA 
and  ORA,  respectively  -  rad 


Not  all  terms  are  necessarily  required;  the  ones  chosen  depend 
on  the  gyro  design  and  on  the  application.  For  example,  some  gyro  engi- 
neers  drop  tho  Doo  term  since  no  satisfactory  physical  explanation  for  it 
has  been  presented,  nor  is  there  incontrovertible  evidence  that  it  exists. 
However,  the  term  has  been  included  for  completeness. 

In  most  gyro  tests,  and  in  strapdown  applications,  the  command 
rate-to-balance  current  is  measured  and/or  recorded.  It  is  usually  assumed 
that  oj^  =  0,  particularly  for  those  gyros  with  synchronous  motors.  This 

a 

assumption  may  not  always  be  valid.  A  common  inperfection  occurring  in 
gyros  is  wheel  hunting,  which  results  from  an  interaction  of  the  rotor 
inertia  with  the  magnetic  stiffness  of  the  motor.  The  combination  acts  as 
a  lightly-damped  torsional  pendulum,  resulting  in  oscillation  of  the  rotor 
angular  rate  about  the  synchronous  value.  The  effect  is  started  by  wheel 
disturbances  and  may  or  may  not  persist  with  time.  Solving  Eq.  (4-1)  for 

s 

the  command  rate-to-balance,  with  »r  =  0,  we  have: 


Kyi  *  "fDF  +  DI*I  +  D0*0  +  DS*S  +  +  DC/OaO+  DSSAS 


+  DIO*IaO  +  DOSaO*S  +  +  “l  +  co  ®S  “  *o 

■  <JS  -  JI>*S  WI 

+  1  r  (4'2) 

r 


Eq.  (4-2)  will  be  assumed  in  all  test  data  reduction  for  obtain¬ 
ing  the  performance  coefficients  L^,  etc. 
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5.  TUMBLE  TESTS 


5.  1  Introduction 


In  the  usual  tumble  tests,  the  gyro  is  mounted  on  a  turntable 
whose  axis  of  rotation  is  nominally  parallel  to  the  earth's  polar  axis.  As 
the  turntable  is  rotated  about  its  axis,  the  earth's  gravity  field  is  rotated 
relative  to  the  gyro  and  many  or  all  of  the  gyro  coefficients  are  exercised. 
The  command  rate-to-balance  current  is  measured  at  a  scries  of  table 
angles  and  recorded.  From  the  data  thus  obtained  ami  auxiliary  information, 
the  gyro  coefficients  are  determined. 

Tumble  runs  are  performed  by  either  of  two  basic  methods: 

a)  Discrete  rotation  runs  in  which  the  table  is  rotated  to  a 
series  of  table  angles  and  the  average  command  rate-to- 
balance  current  is  measured  at  each  table  angle  after 
equilibrium  is  attained. 

b)  The  table  is  rotated  at  a  slow  uniform  rate  (usually 

t  20  x  earth's  rate)  and  the  command  rate-to-balance 
current  is  recorded  continuously.  In  either  method,  the 
average  command  rate-to-balance  current  is  determined 
at  a  number  of  equally -spaced  table  angles  and  the  gyro 
coefficients  are  determined  by  the  method  of  least  squares. 


5.  <2  Error  Sources 


There  are  many  error  sources  in  a  tumble  test,  both  static  and 
dynamic.  Among  the  most  important  are  the  following: 

a)  misalignment  of  the  table  axis  with  respect  to  the  earth's 
polar  axis, 

b)  misalignment  of  the  gyro  reference  axes  with  respect  to  the 
table  axis  and  its  zero  reference  position, 

c)  non -orthogonality  of  the  gyro  output  axis  and  the  spin  axis, 

d)  convection  torques  due  to  density  gradients  (may  be  caused 
by  temperature  gradients  or  by  non-homogeneous  fluid), 

e)  table  axis  of  rotation  changes  with  direction  of  rotation, 
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f)  wobble  of  table  axis, 

g)  torque  signal  lags  due  to  uncompensated  filters  in  the 
torquer  loop  and/or  torque  recorder* 

h)  motion  of  the  gimbal  relative  to  the  gyro  case  due  to  off- 
neutral  buoyancy  temperature  or  unbalances, 

i)  magnetic  fields, 

j)  motion  about  non-principal  axes  of  the  gimbal, 

k)  bubbles,  lint,  or  other  contaminants  in  the  fluid, 

l)  power  supply  variations, 

m)  pier  motion 

n)  noisy  bearings, 

o)  pivot  friction  torque  due  to  off-neutral  buoyancy  tempera¬ 
ture  coupled  with  acceleration. 

We  will  investigate,  in  some  detail,  error  sources  a),  b),  c) 
and  d).  Error  sources  e)  and  f)  can  be  measured,  and  either  corrected  or 
the  data  compensated,  though  either  may  be  difficult  to  accomplish.  Sources 
g)  and  h)  are  dynamic  type  errors  whose  magnitude  and  phase  are  functions 
of  table  rate.  Stray  magnetic  fields  around  the  table  should  be  eliminated, 
as  far  as  practical,  and  the  gyro  protected  by  a  continuous  magnetic  shield 
in  order  to  control  error  source  i).  Error  source  j)  may  be  minimized  by 
careful  gyro  design  and  by  reasonable  care  in  alignment,  as  will  be  shown. 

The  last  five  error  sources  k)  -  o)  generally  cause  discontin¬ 
uous  and/or  erratic  command  rate-to-balance  traces.  Usually,  they  are 
random  in  nature  and  result  in  run-to-run  variations  which  are  frequently 
difficult  to  diagnose.  Bubbles  or  lint  in  the  fluid  may,  under  some 
circumstances,  appear  to  be  a  large,  repeatable  unbalance,  or,  under 
some  circumstances,  they  may  act  like  elastic  restraints,  or  they 
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may  cause  dramatic  shifts  in  the  command  rate-to-balance  measurements. 
Pier  tilts  may  be  eliminated  by  automatic  or  manual  leveling  or  the  tilt 
may  be  monitcred  and  the  data  corrected  (which  may  be  difficult). 

5.  3  Transformation  of  Coordinates 

Let  us  derive  the  general  transformation  equations  required  to 
go  from  inputs  in  an  earth-fixed  coordinate  system  to  a  gyro-case-fixed 
coordinate  system,  including  any  additional  inputs  of  the  table  relative  to 
the  earth.  In  Fig.  5-1,  let  xyz  be  the  earth-fixed  coordinate  system  and 
let  uvw  be  the  coordinate  system  fixed  to  the  gyro  case.  Let  TA  be  the 
table  axis,  which  is  nominally  coincident  with  the  z  axis,  with  the  table 
angle  0  measured  from  the  xz  plane. 


Fig.  5-1.  Earth  and  Gyro  Coordinate  Systems 

The  table  axis  is  slightly  misaligned  with  the  z  axis,  as  defined 
by  the  angles  6,^  and  or.  The  w  axis  of  the  gyro  is  slightly  misaligned  with 
respect  to  TA,  as  defined  by  the  angles  dw^,  and  (3.  Fig.  5- 1(a)  is 
a  view  looking  in  the  minus  z  direction  and  Fig.  5- 1(b)  is  a  view 
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normal  to  the  plane  defined  by  TA  and  z  axeu  with  the  w  axis  rotated  about 
TA  so  that  it  is  in  the  plane  of  TA  and  z.  The  angles  6,pz  and  6^,^  are 
each  «  1  radian.  The  angles  a  and  {3  may  have  any  values  from  0  to  2n. 
Physically,  the  actual  table  axis  may  be  skewed  to  the  z  axis,  then  TA 
in  Fig.  (5-1)  is  taken  to  be  the  projection  of  the  table  axis  on  a  plane  con¬ 
taining  the  z  axis  and  parallel  to  the  actual  table  axis.  Similarly,  the  w 
axis  shown  in  Fig.  (5-1)  is  the  projection  of  the  actual  w  axis  on  a  plane 
containing  TA  and  parallel  to  the  actual  w  axis. 

Lot  x,  y,  /.  be  the  vector  components  of  earth  rate  o"  of 
acceleration  along  the  x,  y  ami  /.  axes,  respectively.  Let  t  be  the  vector 
representing  the  table  rate  relative  to  earth.  Let  u,  v,  ant!  w  be  the  total 
resultant  vector  components  along  the  u,  v,  and  w  axes,  respectively.  The 
vector  components  in  the  uvw  coordinate  system,  due  to  the  inputs  x,  y,  /. 
andT,  are  determined  by  the  following  transformations: 

a)  Rotate  the  xyz  system  about  the  z  axis  through  the  angle 
a  to  form  the  Xj  y^  z  ^  system. 

b)  Rotate  the  x^  y  ^  z^  system  about  the  y^  axis  through  the 
angle  to  form  the  x^  y ^  z %  system.  The  z^  axis  now 
coincides  with  TA,  and  the  vector  t  is  now  introduced. 

c)  Rotate  the  x^,  y  .  v. ^  system  about  the  z,  axis  through  the 
angle  (p  -  u  f  G)  to  form  the  x  ^  y^  ^  system. 

<l)  Rotate  the  x^  y^  /  ^  system  about  the  y  ^  axis  through  tin: 
angle  6w-j>  to  form  the  y^  / ^  system.  The  axis  now 
coincides  with  the  w  axis. 

e)  Lastly,  rotate  the  y^  z^  system  about  the  axis 
through  the  angle  -Pto  coincide  with  the  uvw  system. 

Since  6_  and  6  _  are  <*•  1  radian,  we  may  use  small  angle 
Tz  wT 

approxim ations  for  the  trigonometric  functions,  i.e., 


CO.  6Tz~co.  6wT-  1 


si"  6Tz-«  iTz  and  .in  6  ^  6^ 

6™  •  6  ^  ^  0 
Tz  wT 
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(5-1) 


After  performing  the  above  transformations  and  applying  the 
approximations  of  Eqs.  (5-1),  we  find 

u  x  cos  0  +  y  sin  0  -  1  [ 6  Tz  cos(0  -  *  *  +  6wT  cos  0  ] 

-■*  ®wT  cos  P 

v  *•  -x  sin  0  +  y  cos  0  +  T  [flTz  sin(0  -  «)  -  6wT  sin  P  ] 

-T  fiwT  sinP  (5-2) 

w~  x  [6Tz  cos  or  +  6wT  cos(p  +  0)]  +  y  [6Tz  sin  a 
+  6wT  sin(?  +  0)]  +  z  +7 

5.  4  Continuous  Rotation  Tumble  Test 

We  will  first  consider  the  tumble  test  in  which  the  table  i6 
rotated  continuously  and  uniformly  with  respect  to  earth.  The  discrete 
position  tumble  test  will  be  investigated  in  a  later  section. 

5.  4.  1  Position  No.  1,  ORA  — >  North 

With  the  table  axis  nominally  parallel  to  the  earth's  polar  axis 
EA,  mount  the  gyro  on  the  table  so  that  ORA  is  nominally  parallel  to  TA 
and  points  north*  and  IRA  points  west  when  the  table  angle  0  =  0°.  Fig. 

5-2  shows  the  relative  directior  s  of  the  various  axes.  Fig.  5-2(a)  is  a 
view  looking  south  along  EA  and  Fig.  5-2(b)  is  a  view  normal  to  the  plane 
of  EA  and  TA  with  ORA  rotated  about  TA  into  the  plane  of  EA  and  TA. 

The  misalignment  angles  and  GqnT  correspond  to  6^,z  and  of 
Fig.  5-1  and  are  greatly  exaggerated. 

The  table  will  be  rotated  at  a  rate  boig,  relative  to  earth,  in  a 
counterclockwise  direction  (same  direction  as  earth  rate)  and  then  in  a 
clockwise  direction,  i.  e. ,  -btig,  relative  to  earth. 
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Fig.  5-2.  Position  No.  1 

5.4. 1.  1  Counterclockwise  Rotation  of  the  Table;  The  angular  velocity 
components  to  be  substituted  in  Eqs.  (5-2)  are  x  *  0,  y  =  0,  z  =  and 
7  a  b«£.  The  angular  velocity  vectors  b^,  Bq(  and  es  correspond  to  the 
vectors  u,  w,  and  -v  of  Eqs.  (5-2),  respectively. 

Bj  -  -b£  [(1+b)  OqnT  CO#  ^1  +  ®  TE  co,(0-*)3 

•o“»E<1+b> 

®S  “e  t(1+b>  fiONT  ,ln  *1  ■  fiTE  •in(6-«)]  (5-3) 

•o‘° 

The  acceleration  components  to  be  substituted  in  Eq.  (5-2) 
are  x  =T  *0,  y  *  cosX,  z  =  sinX,  where  X  is  the  astronomic  latitude  (positive 
in  northern  hemisphere).  The  acceleration  components  a^,  Sq,  ag  correspond 
to  the  vectors  u,  w,  and  -v,  respectively  of  Eqs.  (5-3). 

aj  »«  cosX  sine  -  sinX  WqnT  co#  *1  +  °TE  co*(0'®)3 
»q  =*■  sinX  +  cosX  C&ONT  ,in^l  +  0)  +  ®te  ,lnflf  3 

ag  “■  -cosX  cos©  +  sinX  C80NT  ,ln  ^1  "  ®TE  •ln(0'o,)3 
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At  middle  latitudes  (25°  to  60°)  where  practically  all  gyro 
testing  is  performed,  the  effects  of  small  misalignments  on  the  accelera¬ 
tion  components  will  result  in  only  second-order  errors  in  the  drift  rates 
and  will  be  neglected. 

a^  °»  cosX  sin0 

aQ  =*•  sinX  (5-4) 

ag'S'-cosX  cos© 

Substitute  Eqs.  (5-3)  and  (5-4)  in  Eq.  (4-2),  Let  Iq^  be  the 

command  rate -to- balance  current  in  Position  No.  1  (ORA  — >  North)  with 

counterclockwise  rotation  at  table  angle  0.  =  k0  ,  where  0  =  360°/n 

k  n  n 

and  k  =  1,  2,  3, . . . ,  n.  The  terms  are  arranged  as  in  a  Fourier  series. 

KT*ONk  =  +  DC  8in*  +  D00  8in?X  +  1/2(Dji  +  cos2x 

+  u»E  (i  +  b)  6qnt  (cospj  -  :Q  sin  Pj) 

JS  '  JI  2  , ,  .  1_v2  ,2 
412,600  Hr  *E  (1  b)  6ONT  ln2PlJ 

-  CDj  cosX  +  1/2  Djq  sin  2X  +  (»ino+  GqCOS 

*  “’E  (1  +b)  flONT  6TE  cos(o+P1)]sinken 

+  l/2[Dlscos2x-Zg|>30-oh  4  62Ecos20r  jsin2ken 

r 

+  [Dg  cosX  +  1/2  DQg  sin2X-  6TE  (cosor  -  Cq  sinar) 

•20M60H1  4  (1  +  b)  6ONT 6TE  8i“(®+Pi)]cosken 
*  r 

+  1/2[(DU  -  Dgg)  cos2X-  4  6TEsin2°^COS2k 

f  (5-5) 
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5.4. 1.2  Clockwise  Rotation  of  the  Table;  The  setup  is  the  same  as  in 
Fig.  5-2.  but  the  table  is  now  rotated  in  the  opposite  direction  without  dis¬ 
turbing  the  setup  used  for  the  counterclockwise  rotation.  The  command 
rate-to-balance  equation  is  the  same  as  Eq.  (5-5).  except  that  (1  +  b)  is 
replaced  by  (1  ~  b).  and  the  command  rate-to-balance  current  is  1q^  at 
table  angle  6^  *  k$n.  Though  the  table  is  rotating  clockwise,  the  table  angle 

6,  is  still  measured  in  a  counterclockwise  direction  from  the  west,  as 
k 

shown  in  Fig.  5-2(a). 

5. 4. 1.  3  Weighting  the  Data  from  Position  No.  1:  In  order  to  eliminate 
some  of  the  effects  of  the  gyro  misalignment  angle  and  to  minimize 

the  remaining  effects,  we  will  weight  the  data  as  follows: 

WONk  *  "ZTj"  lONk  +  (b  +  1)  lONk^  *  "^DF  +  D0  8inX 

+  Dqq  sin2X  +  1/2(DU  +  Dgg)  cos2X 

+  412%00ir  “e  6ONT  *ln  2pl^’  ^DI  co,X 
+  1/2  DJO  8in2X  +  »E6TE  (sin®  +  C0  cos®)]  sin  kGn 

+  1/2  [Djg  cos2X  +  20fe,  TWIT  *4  °TE  CO®  20r>in  2ken 
+  [Dg  cosX  +  1/2  Dgg  sin  2X  -  «»E  6TE  (cos® 

-  CQ  sin®)]  cos  k6n  +  / 2  [(Dn  -  Dgg)  cos2X 

’  2W,  W  Hr  *E  ®TE  8ln  2ar]  008  2kGn  (5'6) 

5. 4. 2  Position  No,  2,  ORA  — >  South 

Position  No.  2  is  obtained  from  Position  No.  I  by  rotating  the 
gyro  180°  about  SRA  so  that  ORA  points  south  along  TA  and  IRA  points  east 
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when  9  =  0°.  Rotating  about  SRA  minimizes  the  time  required  for  settling. 

In  the  following  analysis,  it  is  assumed  that  the  position  of  TA  relative  to  EA 
is  unchanged.  If  this  is  not  true  in  any  particular  test  setup,  then  the  analy- 
sis  must  be  appropriately  modified  and  becomes  more  complicated.  However, 
in  general,  the  angle  that  ORA  makes  with  TA  would  not  be  the  same  as  in 
Position  No.  1. 

Fig.  5-3  shows  the  relative  directions  of  the  axes.  Fig.  5-3(a) 
is  a  view  looking  south  along  EA,  and  Fig.  5- 3(b)  is  a  view  normal  to  the 
plane  of  EA  and  TA  with  ORA  rotated  about  TA  into  the  plane  of  EA  and  TA. 


Fig.  5-3.  Position  No.  2 


5.  4.  2.  1  Counterclockwise  Rotation  of  the  Table:  Let  the  table  rotate 
counterclockwise  (0  >  0)  at  a  constant  rate  bu>£  with  respect  to  earth.  The 
angular  velocity  vectors  to  be  substituted  in  Eqs.  (5-2)  are  x  -  “  =  0,  “  =  n 
and  t  =  buij,.  The  angular  velocity  components  o>  ,  and  u)r  correspond 

to  -u,  -w,  and  -v,  respectively  of  Eqs.  (5-2). 
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•j  *•  -«»E  C  (1+b)  6osT  co.  p2  -  6te  co.(e-or)] 
wo"mU,E  (1+b) 

r  (5-7) 

"S  “■  WE  ^1+b)  ^OST  ,in  *2  “  ®TE  *in 

®o*‘  0 

The  equivalent  applied  acceleration  component,  along  the  gyro 
reference  axes  are  also  determined  from  Eq%  (5-2)  and  simplified. 


*■  -co.X  .in0 


aQ  *■  -sin\ 


(5-8) 


-co.X  co.6 


Substitute  Eqa.  (5-7)  and  (5-8)  in  Eq.  (4-2).  Let  i^g^  be  the 
command  rate -to- balance  current  in  Position  No.  2  (ORA  — >  South)  with 
counterclockwise  rotation  at  table  angle  0^  =  k6R»  where  k  =  1,  2,  3, ....  n. 

KTiosk  =  -£°f  "  Do  8inX  +  Doo  "k*2*  +  */2(Dn  +  Dss^  co,2x 

+  wE  (1+b)  fiQST  (co.  P2  -  CQ  «in  P2) 

+  *i77Wrr  “e  (1+b)2  °OST  Sin  2P2] 

+  [Dj  co.X  -  1/2  Djq  .in  2X  +  '"e^TE  (8^n®  "  Coc°8*) 

+  206,  W  H  ®E  (l+b)  6OST  6TE  c°8<*+P2^  8in  k0n 

-  1/2  [Djg  cos2X  +  206,  WlT  “e  6TE  c°8  8ln  2k6n 
+  [Dg  co.X  -  1/2  Dqq  sin  2X  +  u>E  6TE(cosor  + CQ  sincr) 
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+  20?,  W Hr  4  (1+b)  flOST  6te  •i“(*-P2)]co8  ken 

+  1/2[(DU-DSS)  co.2X  +  Zoi,mWOU  4  4E8in20,^CO82ken 

r  (5-9) 


5. 4.2.2  Clockwise  Rotation  of  the  Table:  The  setup  is  the  same  as  in  Fig. 
5-3,  but  the  table  is  row  rotated  in  the  clockwise  direction  without  disturbing 
the  setup  used  for  the  counterclockwise  rotation  of  the  table.  The  command 
rate-to-balance  equation  is  the  same  as  Eq.  (5-9),  except  that  (1+b)  is  re¬ 
placed  by  (1-b),  and  the  command  rate-to-balance  current  is  i^g^  at  table 
angle  0k  =  kG  .  The  table  angle  0  is  still  measured  in  a  counterclockwise 
direction  from  west  as  shown  in  Fig.  5-3(a). 


5.  4.  2.  3  Weighting  the  Data  from  Position  No.  2:  In  order  to  eliminate 
some  of  the  effects  of  the  gyro  misalignment  angle  an(*  to  minimize 

other  misalignment  effects,  we  will  weight  the  data  as  follows: 


K„ 


W 


OSk  =  TF  £(b-1)  lOSk  +  (b+1)  iOSk^ 


=  -[Df-Dq  sinX+D00  sin2\  +  l/2(Du+DgS)  cos2X 


Jc.  -  J 


+  TI775MTT  4  (fe2"l)  6ost  8in  2p23  +  tDi  cosX 


-  1/2  Djq  sin  2X  +  u>E  Gaisina  -  Cq  cosor)]sin  kGn 

-  1/2  [Djg  cos2X  +  206^  300  Hr  *E  6TE  COS  Za  ^  8in  2kGn 
+  [Dg  cosX  -  1/2  DQg  sin  2X 

+  UJj,  (cosor  +  CQ  sin«)]cos  +  l/2[(DIj-Dgg)  cos2X 

*  mfwrr;  4  6TE  3in  20f  1  008  2ke  (5“lu) 
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5.4.3 


Control  of  Misalignment  Errors 


It  is  somewhat  easier  to  investigate  the  impact  of  the  misalign¬ 
ment  errors  if  we  look  at  the  sum  and  difference  of  Eqs.  (5-6)  and  (5-10). 

WONk  +  WOSk  =  *^2DF'f  2DOO  8i“2x  +  fDII+DSS*  co®2* 

+  TOTTO'Hr  *E  (l>2*1)  (6ONT  8in  2pl  +  6OST8ln  2p2n 

-  [DJO  sin  2X  +  2  «E  6TE  Cq  cos®]  sin  k0n 

+  2[Dg  cosX  +  (Bg  6te  Cq  •in®]  cos  k0R 

+  (Djj  -  Dgg)  cos2X  cos  2k0n  (5-11) 

WONk  *  WOSk  *  -[2  D0  +  4lTrmnt  4tb2-l)(#|NT.in2(l 

-  ®0ST  *in  2P2)]  -2[Dj  cosX  +  u^,  6TE  sin®  ]sin  k6n+p>IS  cos2X 

♦  206!  M0ffr  4  6TE  008  2or]  8in 2ken  +  ^OS  8in  2X 

-  2we6te  co80flc°8  k0n  *  zsItWh;  4  6te  8ln  20f  co#  2ken 

(5-12) 

It  is  easier,  in  general,  to  control  the  misalignment  of  the  table 
axis  to  the  earth  axis  (angle  O^g)  than  it  is  to  control  the  misalignment  of  the 
gyro  output  reference  axis  with  respect  to  the  table  axis  (angles  Oqj^j.  and 
Oqst^*  The  latter  misalignment  angles  result  from  manufacturing  and 
assembly  tolerances  of  the  gyro,  the  mounting  fixture  and  the  table.  In 
addition,  those  angles  may  be  a  function  of  the  clamping  loads,  the  average 
temperature,  and  the  temperature  gradients  in  the  test  setup  and  in  the  gyro. 
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It  is  relatively  easy  to  control  0Tg  to  less  than  0.  001  rad  (“*3  min),  and  it 
may  be  controlled  as  closely  as  15  x  10"^  rad  (**3  sec)  on  a  very  stable 
pier. 

The  misalignment  angle  Cq  in  *  well-designed  and  assembled 
gyro  should  be  less  than  0.01  rad  (**34  min);  therefore,  if  we  control  u>E 
in  Eq.  (5-12)  to  a  satisfactory  degree,  then  the  error  terms  containing 
u)g  flrpg  Cq  In  Eq.  (5-11)  should  be  quite  negligible  and  may  be  deleted. 

At  a  latitude  of  45°,  it  is  seen  that  tUg  Oje  mu8t  be  l*88  than  7% 
of  Dj  and  5%  of  Dq<.  if  these  coefficients  are  to  be  determined  with  an  error 
of  less  than  10%.  For  example,  if  Dq^  =  0.  001  nig,  then  2  uig  0^,g  must  be 
less  than  0.  1  x  0.  001  o>g  or  0^,g  must  be  less  than  1  x  10"*/2  =  5  x  10*^  rad 
(*■  10  sec  >ifDos  is  to  be  determined  with  an  error  of  less  than  10%.  It  is 
obviously  far  more  difficult  to  determine  Dq<<  or  Dj  accurately  if  they  are 
only  1  x  10  u>g,  though  some  manufacturers  make  even  more  doubtful 
claims. 

Let  us  now  look  at  the  remaining  error  terms;  these  involve  the 
ratio  (Jg  -  Jj)/Hf.  Since  iu  general,  the  misalignment  angles  6qNT  and 
6qst  would  be  greater  than  6^,g  and  (b^-l)  » 1,  then  the  errors  in  the  con¬ 
stant  terms  of  the  Fourier  expansions  are  more  significant  than  the  errors 
in  the  other  Fourier  coefficients.  The  Central  inertial  Guidance  Test 
Facility  at  Holloman  AFB  normally  performs  its  tests  with  b  *  20  and,  on 
occasion,  as  high  as  b  =  100.  Let  us  look  at  the  order  of  magnitude  of  these 
errors  by  assuming  some  valuer  for  the  various  parameters  which  will  give 

errors  greater  than  what  would  ordinarily  be  expected.  Let  Je  -  =  1000 

2  4  2  o  i 

gm*  cm  ,  Hr  =  10  gm*  cm  /s,  6qNT  *  Oqsx  =  0,05  rad»  sin  2^1  =  sln  2P 2 

=  0.  5,  b  =  100,  and  o>g  15  deg/hr.  Then, 

?T^Wir  4  (b2‘**  (4ONT  *ln  2(1  *  *OST  ,ln  = 

412,608x10. 000  x  225  x  l0- 000  0025  x  °- 5 

+  0.  0025  x  0.  5)  “■  0.  0014  deg /hr 
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Therefore,  even  with  these  rather  extreme  values,  the  error  would  be 
negligible  in  most  applications.  However,  it  would  be  well  to  estimate  this 
error  for  each  new  design  and  test  setup  in  order  to  confirm  that  it  is 
negligible.  If  the  error  is  not  negligible,  there  are  several  steps  we  may 
take  to  reduce  it;  they  are:  a)  apply  closer  tolerances  on  the  manufacture 
and  assembly  of  the  gyro,  b)  apply  closer  tolerances  on  the  manufacture 
and  assembly  of  the  mounting  fixture,  c)  correct  any  non-orthogonality  of 
the  table  face  to  the  table  axes,  d)  use  a  lower  table  rate,  e)  increase  the 
angular  momentum  of  the  gyro  and/or  decrease  Jg  -  Jj. 

5.  4.  4  Gyro  Coefficients  from  OA  — >  North  and  OA-»  South  Tests 

If  we  assume  that  the  misalignment  terms  of  Eqs,  (5-11)  and 
(5-12)  are  controlled  and  negligible  compared  to  the  gyro  coefficients,  as 
is  necessary,  and  we  now  apply  the  method  of  least  squares  to  Eqs.  (5-11) 
and  (5-12),  then  the  following  equations  are  obtained: 

n 

2  DF  +  2  Dqo  sin  X  +  <Du  +  Dss*  C°*  ^^WONk+  WOSk^ 

DIO  =  n  slI'2X  (WONk  +  WOSk)  ,ln  ^n 

DS  =  n  cosX  k^i  (W°Nk  +  '^OSk*  co*  k0n 

DII  "  DSS  =  ^ONk  +  WOSk*  008  2k0n 

n  cos  A  K=1 

D0  =  2n  sinX  ^  (WONk  ’  WOSk) 

-1  n 

D.  - - —r  T  (WAXI,  -  W^J  sin  k0 

I  n  cosX  '  ONk  OSk'  n 
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(5-13) 


=  7^77  I,  ^  "» 

DQS  3  ^  (WONk  "  WOSk)  coa  k0n 

Since  there  ere  only  two  equations  Involving  the  four  unknown 
coefficients  Dp,  D^,  DqO’  end  Dgg,  it  is  necessary  to  obtain  two  more  inde¬ 
pendent  equations  involving  Dp  and  at  least  one  of  the  other  three  coefficients, 
or  make  some  simplifying  assumptions. 

Seme  organizations  assume  that  Dqq  3  0,  in  which  case  only  one 
more  independent  equation  is  required.  An  additional  tumble  test  with  SRA 
parallel  to  the  earth's  polar  axis  will  provide  the  necessary  relationship. 

The  equations  for  such  a  tumble  test  are  derived  in  a  manner  similar  to 
that  used  for  Positions  Nos.  1  and  2.  Unfortunately,  a  new  misalignment 
error  term  appears  which  involves  Jq/H These  errors  are  probably 
negligible  in  most  cases,  but  if  not,  they  are  not  so  readily  amenable  to 
minimization  by  design  changes. 

Other  organizations  assume  that  Dqq  and  either  D^j  or  Dgg  are 
zero.  Equations  (5-13)  are  then  sufficient  to  determine  all  the  remaining 
coefficients.  The  assumption  that  Dqq  *  0  is  justifiable  on  the  somewhat 
weak  grounds  that  no  theoretical  analyses  have  shown  a  need  for  such  a  term. 
However,  it  is  difficult  to  justify  setting  either  Djj  or  Dgg  equal  to  zero, 
and  certainly  not  one  in  preference  to  the  other. 

The  most  accurate  method  for  obtaining  the  quadratic  and  cross- 
coupled  coefficients,  such  as  D^  and  Dgg,  is  on  a  precision  linear  shaker. 
The  common  electromagnetic  shakers,  used  either  with  or  without  a  slip 
table,  are  not  usually  satisfactory  since  they  may  have  substantial  cross 
axis  and  angular  motions.  This  technique  will  be  analyzed  in  Section  6. 

We  will  now  turn  to  the  discrete  position  tumble  test. 
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Discrete  Position  Tumble  Test 


S.  5 

In  the  discrete  position  tumble  test,  the  command  rate-to- 
br.iance  current  is  measured  at  a  series  of  table  positions  after  equilibrium 
has  been  reached*  Thus,  some  of  the  dynamic  error  terms  are  avoided,  but 
of  course,  other  errors  appear. 

5.  5. 1  Position  No*  1*  ORA  — >  North 

The  mounting  position  is  identical  to  that  described  in  Section 

5. 4.  1  for  the  continuous  rotation  tumble  test.  The  command  rate-to-balance 

current  is  measured  at  table  angles  0,  =  k0  where  0  =  360°/n  and 

k  n  n 

k  =  1,2,  3 . n.  The  results  are  the  same  as  given  by  Eq.  (5*5)  with  b  =  0. 

The  gyro  must  be  allowed  to  attain  equilibrium  at  each  table  angle. 

^T^ONk  s  “^DF  +do  8*n*  +  ^OO  8*°  *+  1/2<Dn+DSS*  008  * 

+  ®ONT  ^co8*l  *  8ia  ^ 

Jg  ”  ^1  2  2 

+  417, '600  Hr  “e  flONT  8ln  2pl^ 

-  [Dj  cosX  +  1/2  DIO  sin  2X  +  (Og.  0TE  (sincr  +  CqCOSot) 

’  206,  300  Hr  ®E  °ONT  ®TE  cos<*+  Pf)]  ,in  ken 

+  1/2  [Djg  cos2X  +  2M'  <4  flTE  c°8  ZeJBtoZkQn 

+  [DscosX+  1/2  Dqq  sin  2X-  <ke6tE  (coso-  Cq  sino) 

'  ZoI/WfT  4  °ONT  6TE  Sin(a+  Pl)]  COS  k0n 
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(5-14) 


+  1/2  [Dn  -  Dgg)  cos2X 

•  •*  6te  8in  2of]  008  2k^ 

Unlike  the  continuous  tumble  teste,  it  makes  no  difference, 
theoretically,  whether  the  table  is  rotated  clockwise  or  counterclockwise 
from  one  table  angle  to  the  next,  providing  equilibrium  conditions  are 
attained  at  each  table  angle  before  readings  are  taken,  and  providing  there 
is  no  hysteresis  due  to  pivot  friction  or  other  causes.  If  hysteresis  effects 
are  suspected,  then  readings  should  be  taken  in  both  clockwise  and  counter¬ 
clockwise  rotation.  The  rotation  from  one  position  to  another  should  be 
smooth  and  unidirectional. 

5.5.2  Position  No.  2,  ORA  — >  South 

The  mounting  position  is  identical  to  that  described  in  Section 

5.  4.  2.  The  command  rate-to-balance  current  is  measured  at  table  angles 

0,  =  k0  t  where  9  =  360°/n  and  k  *  1,2,  3,,..,  n.  The  results  are  the 
k  n  n 

same  as  given  by  Eq.  (5-9)  with  b  =  0.  The  same  cautions  regarding 
equilibrium  and  hysteresis  given  in  Section  5. 5. 1  also  apply  here. 

KTlOSk  *  "tDF  -  DQ  sinX  +  Dqq  sin2X  +  l^Djj+Dgg)  cos2X 
+  *E  °OST  (c°8  h  "  co  8ln  P2) 

+  412/600  Hr  4  6ost  8in  2p2]  +  [DI  cosX 

-  1/2  DIO  sin  2X  +  u>E  0TE  (sinar-  Cq  cosar) 

+  Z06/300H  ®E  °OST  6TE  co,(*+  8ln  ken 
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-  1/2 [Djg  cos2X  +  “e  6te 008  2*^ 8ln  2k6n 

+  [Dg  co*X  -  1/2  Dos  sin  2X  +  wE  0TE  (cos*  +  Cq  sin*) 

+  Zfl5,nroHr  ®E  *OST  ®TE  8in*®  "  p20c°®  ken 

+  1/2[(Djj-  Dgg)  cos  X  +  25^7300  FT  Wg  6TEsin2»]cos  2k0n 

r  (5-15) 

5.  5.  3  Gyro  Coefficients  from  OA  — >  North  and  OA  — >  South  Tests 

The  sum  and  difference  of  Eqs.  (5-14)  and  (5-15)  are: 

*Sr^ONk  +  *osk^  *  -£2Df  +  2Doo  8ln  *  +  ^DII+DSS^  COS  X 
+  ®e  ®ONT  ^CO*  ^1  “  co8ln 
+  WE  ®OST  ^CO*  ^2  "  ^ 

+  4TZ7WTTH  *E  *ftONT  Bin  2p  1  +  ftOST  8in  2  P2^ 
r 

-  |Diq  sin  2X  +  2^  6te  Cq  cos  * 

“  1557351? Hr  “e  °te  ^6ont  co®(®+  Pi> 

+  6 QgT  cos(*  +  P2)  ]  |  Sin  ken 
+  |2DS  cosx  +  2  *E  6te  Cq  sin* 

'I557555Tr  ®E  6TE  C6ONT  8ln  *1* 

-  6qST  #ln<®“P2)  3 1  cos  R0q  +  ^DU”  DSS^co*2*  008  2k6n 

(5— 1.6) 
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^^ONk  "  *OSk^  3  ’  t2D0  ,ln^  +  ®E  ®ONT  ^CO*  ^1*^0  ,in  M 


“  *E  °OST  *COi  *2  "  ^Oiln  irySOOTT,'  "fc^ONT  ,ln2pl 
-62)ST«ln  2|2)]  -  |2  Dj  co«X  +  2  «E  6TE  »lnflf 

Js-Ji  2 

*  m;mur  we  6te  [fioNT  co,(*+ 

-  6q§t  coi  (or  ♦  l2)]|  sin  k6n 

+  [Djg  co«2X  +  g fly,  “E  ®TE  COS  2»3®in  2k0n 

+  |^OS  2*  “  2  ®E  ®TE  co,0f 
Jg  "  2 

"  ZW,  300’Hr  “E  6TE  t6ONT  ,in(0,+  pi) 

+  6OST  •i“(*-P2)3|co*  ken 

■  TB^WTT  *E  ®TE  ,ln  2*  CO*  2ken  (5"17> 

The  discussion  in  Section  5.  4.  3  on  the  error  terms  containing 

(Jc  -  Jt)/H  would  also  apply  here,  except  ihat  the  errors  are  usually 
sir  £ 

negligible  since  they  are  not  multiplied  by  (b  -1).  It  is  obvious  that  any 
remaining  error  terms  containing  products  of  misalignment  angles  are  very 
easily  controlled  and  may  also  be  neglected.  Neglecting  those  terms, 

Eqs.  (5-16)  and  (5-17)  reduce  to: 

^T^ONk  +  iOSk^  3  "  ^2Df  +  2Doo  ,in  *  +  ^DII  +  DSs3  cos  X 

+  ®E  ^®ONT  COB  ^1  +  ®OST  co*  ^2^  "  DIO  *in  2*  sin  k®n 
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+  2Dg  cosX  cos  k0n  +  (Djj  -  Dgg)  co*2*  COB  2k0n  (5-18) 

^T^'  ONk  “  lOSk^  *  -E2D0  ilnX  +  “E^ONT00*  *1“  ®OSTCO*  *2^ 

> 

•  2(D|  coaX  +  »E  6^e  ainor)  a  in  k0R  ♦  D-gCoa  X  Bin  2kGn 

+  (Dqs  ain  2X  -  2»E  flTE  coacr)  cos  k0n  (5-19) 

Knowing  the  deaired  accuracy  for  the  performance  coefficients, 
it  is  relatively  easy  to  control  the  misalignment  angle  6^.^  to  the  required 
magnitude.  Unfortunately,  the  terms  ci>E  co*  ^  end  »E  6 qst  co*  ^2 

are  not  as  amenable  to  control  and  cannot  be  eliminated  by  weighting  clock¬ 
wise  and  counterclockwise  rotation  as  is  done  in  continuous  rotation  tumble 
tests.  A  careful  estimate  of  the  maximum  error  must  be  made  and  compared 
to  the  allowable  error.  The  misalignment  angles  include  misalignment  of  the 
gyro  to  the  table  axis,  non-orthogonality  of  table  face  to  table  axis,  mount¬ 
ing  fixture  errors,  effects  of  bolting  torques,  and  dirt  on  mounting  surfaces. 

Assuming  the  misalignment  errors  have  been  properly  controlled, 
we  may  apply  the  method  of  least  squares  to  Eqs.  (5-18)  and  (5-19).  The 
resulting  normal  equations  are: 

-  K  n 

2DF  +  2D00  *ln2x  +  *Dn  +  DSS)co#  X  *  IT  ^  ^iONk+iOSk) 

DIO  *in  2X  3  jEj  ^ONk  +  4OSk^  8in  k0n 

DS  co,X  *  TT  X,  <lONk  +  lOSk>  co*  k8n 
k*l 

*DII  "  DSS*  CO *  X  s  “n”  jEj  ^lONk  +  lOSk*  COS  2k8n 
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D0  ,lnX  *  “ST  X,  *lONk  *  *08^ 


k*l 


DI  co>X  =  “IT  Xj  (4ONk  "  lOSk}  8in  k0n 


2  2KT  n 

DIS  C08  X  =  “n"  Xj  <lONk  '  *OSk^  ,ln  2k0r 


2K_  n 

DOS  8ln  2X  *  “n~  ,X.  ^ONk  "  iOSk)  cos  k0n 

k-1 


As  with  the  continuous  rotation  tumble  test,  there  are  ten  unknown 
coefficients,  but  only  eight  equations.  This  leads  to  the  need  for  two  more 
independent  equations  or  for  some  simplifying  assumptions.  Refer  to  Section 
5.4.4  for  a  discussion  of  this  subject. 
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6.  PRECISION  LINEAR  VIBRATION  TEST 


6. 1  Introduction 

As  previously  mentioned,  the  moat  accurate  method  of  obtaining 
the  coefficients  of  the  quadratic  and  cross -coupling  terms  such  as  and 
D^g,  is  on  a  precision  linear  vibration  table.  The  great  advantage  of  the 
vibration  test  over  the  tumble  test  is  due  to  the  higher  acceleration  level 
available  which  greatly  magnifies  the  response  because  the  drift  rate  varies 
as  the  square  of  the  input  acceleration.  Since  it  is  the  quasi-static  value 
of  each  coefficient  that  is  generally  desired,  care  must  be  taken  to  have  the 
input  frequency  well  below  any  resonant  frequency  of  the  gyro  or  of  the 
mounting  fixture.  In  research  and  development,  it  is  probably  desirable 
to  test  at  more  than  one  frequency. 

Mounting  the  gyro  and  its  fixture  directly  on  the  ordinary  electro¬ 
magnetic  shaker,  or  on  an  ordinary  slip  table  driven  by  such  a  shaker,  i~ 
often  not  acceptable  because  of  the  large  components  of  cross  acceleralio:. 
and  angular  acceleration  which  ara  usually  present,  A  precision  linea* 
shaker  such  as  the  one  in  the  Charles  Stark  Diaper  Laboratory,  driving  - 
good  Team  table  is  required. 

In  order  to  minimize  the  effects  of  earth's  rate,  the  vibration 
input  should  be  along  a  line  normal  to  the  polar  axis  such  as  a  horizontal 
east-west  line,  and  IA  should  be  normal  to  the  polar  axis  whenever  possible. 
We  will  assume  that  reasonabl  care  has  been  taken  in  the  setup  and  mount¬ 
ing  of  the  gyro  so  that  the  total  misalignments  of  the  gyre  axes,  as  mounted 
with  respect  to  the  vibration  axis,  earth  axis,  and  the  gravity  vector  are 
less  than  0.01  rad  (“■  34  min)  and  may  be  neglected. 

The  determination  of  only  three  coefficients  will  be  illustrated 
in  the  following  sections;  all  others  may  be  determined  in  a  similar 
method. 


Preceding  pita  blink 
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.2 


Determination  of  -  Position  No.  1 


Mount  the  gyro  on  the  precision  linear  shaker  with  ORA  pointing 
north  along  the  earth's  axis  EA,  and  IRA  pointing  east  as  shown  in  Fig.  6-1. 
The  line  H-H  represents  the  horizon,  X  is  the  astronomic  latitude,  and  is 
the  vibration  input  in  £  units  with  amplitude  Aj^  (0  to  peak). 


(a)  VIEW  LOOKING  WEST  (b)  VIEW  LOOKING  SOUTH  ALONG  EA 


Fig.  6-1.  Gyro  Axes  in  Vibration  Position  No.  1 


The  gyro  is  vibrated  at  (n+1)  different  acceleration  amplitudes 
Alk*  •  •  • » *»d  A  jq  =  0.  The  components  of  acceleration  and  earth 

rate  along  the  gyro  axes  in  Position  No.  1  are: 

*m  *  Aik ,ln  « 

•oik  *  ,,bX 


*Slk  *  co,x 


•lE  “  *SE  “  0 
“OE  *  *E 


(4-1) 
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Substitute  Eqs.  (6-1)  and  (6-2)  in  Eq.  (4-2)  and  let  i^  be  the 
command  rate-to-balance  current  in  Position  No.  1  at  vibration  amplitude 

Alk* 

Wjik  *  Kj.  iik  *  -LDp  +  DjAlk  sin  u>t  +  DQ  sinX  +  Dg  cosX 

2  2  2  2 
+  Djj  Ajk  sin  »t  +  Dqq  sin  A  +  Dgg  cos  X 

+  DlQAlk  sinX  sin  at  +  D^g  sinX  cosX  +  Djg  cosX  A^sin  at] 

The  average  command  rate-to-balance  over  a  time  interval  NT. 
where  N  is  an  integer  and  T  =  2 n/a  is  the  period  of  the  vibration,  is: 


W 


K„ 


NT 


Ilk3  W 


f  ‘ikdt*- 

J  A 


Cdf  +  d 


2  2 

+  dOO  ,ln  x  +  Dss  co*  X  +  D 


q  sinX  +  Dg  cosX  +  1/2DQA 
Qg  sinX  cosX]  (6-3) 


2 

lk 


Perform  the  above  test  for  k  »  1,  2f  3, . . . ,  n.  In  testing  a  new 
model  of  a  gyro  or  a  gyro  mounted  on  a  new  fixture,  it  is  good  practice  to 
test  at  several  frequencies  in  order  to  ascertain  that  the  final  test  runs  are 
not  performed  near  or  beyond  any  resonant  frequency  of  gyro  or  mount. 

Without  disturbing  the  setup  in  any  manner,  determine  the  aver¬ 
age  command  rate-to-balance  current  with  Ajk  =  A ^  =  0.  The  average  com¬ 
mand  rate-to-balance  is: 

K  NT 

Who  “ht  f  lio  dt  3  "^df  +  Do  ,lnX  +  Dsco*x  +  Doo  8in  x 

Jo 

+  Dgg  co.2X  +  Dqq  sinX  cosX  ]  (6-4) 

Subtract  Eq.  (6-4)  from  Eq.  (6-3)  for  each  input  acceleration  level, 
sum  the  results,  and  solve  for  Dn- 
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(6-5) 


£ .  Aik  =  2n  "ho 

k*l 


-2  £  W 


k=  1 


Ilk 


It  is  often  useful  to  plot  the  change  in  indicated  command  rate  (Wjjk“wjjQ) 
versus  input  acceleration  level  squared.  This  relationship  should  be 
linear,  and  the  slope  is  proportional  to  D^. 


6.3 


Fig.  6-2. 


Determination  of  Dgg  -  Position  No.  2 

o 

Rotate  the  gyro  90  about  ORA  to  Position  No.  Z,  as  shown  in 


a  =  la 


EA  (North) 


oy  =  A  sinut 


ORA&EA 


(a)  VIEW  LOOKING  WEST  (b)  VIEW  LOOKING  SOUTH  ALONG  EA 


Fig.  6-2.  Gyro  Axes  in  Vibration  Position  No.  2 

Perform  the  vibration  tests  at  n  amplitude  levels  (not  nec¬ 
essarily  the  same  levels  used  in  Position  No.  1,  except  that  A2Q  =  Ajq  =  0). 
The  components  of  acceleration  and  earth  rate  along  the  gyro  axes  are: 


*I2k  *  *C°8X 
a02k  =  ,inX 

aS2k  =  A2k  §ln  •*  (6’6) 


®IE  “■  *SE  “  0 
®OE  -  *E 


(6-7) 
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Substitute  Eqs.  (6-6)  and  (6*7)  in  Eq.  (4*2)  and  reduce  the 
data  in  a  manner  similar  to  the  methods  used  in  Section  6. 2. 

DSS  ^  A2k  *  2n^I20  ‘  2  ^I2k  (6“8) 


6. 4  Determination  of  Djg  •  Position  No.  3 


Rotate  the  gyro  about  ORA  to  the  position  shown  in  Fig.  6*3. 


(a)  VIEW  LOOKING  WEST  (b)  VIEW  LOOKING  SOUTH  ALONG  EA 


Fig.  6-3.  Gyro  Axes  in  Vibration  Position  No.  3 
The  components  of  acceleration  and  of  earth  rate  along  the  gyro 

axes  are: 

\rZ  \Z 

*I3k  =  "  T  C0BX  +~7  A3k  ,in  ®t 


*03k 


=  sinX 


*S3k 


cosX  + 


sin  et 


"iE-  ®SE 
®OE  3  *E 


(6-9) 


(6-10) 
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Again,  take  data  at  each  amplitude  level  A^.  including  A^q  *  0. 
Applying  the  same  data  reduction  procedure  aa  before,  we  find: 

<DU  +  DSS  +  DIS»  i.A3k  *  4“ *130  -  4  X  *I3k  <6-U> 

k=l  k=l 

One  could  substitute  the  values  of  and  Dgg  obtained  from  Eqs. 
(6*5)  and  (6-8)  in  Eq,  (6-11)  and  solve  for  D^g.  Because  of  the  standard 
errors  of  the  mean  values  of  and  Dgg,  it  is  a  somewhat  better  procedure 
to  rotate  the  gyro  90°  about  ORA  from  position  3  and  perform  another  series 
of  vibration  runs  with  the  vibration  inputs  along  IRA  and  SRA  being  180°  out 
of  phase.  The  reduced  equation  for  ttis  fourth  position  is: 

(DII  +  DSS  *  DIS)  A4k  a  4n^I40  *  4  j^j^Mk  (6-12) 

The  value  of  D^g  is  obtained  from  Eqs.  (6-11)  and  (6-12).  The 
coefficients  DIO  and  may  be  determined  in  a  similar  manner. 
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7.  TORQUES  DUE  TO  DENSITY  GRADIENTS 


If  the  gyro  fluid  has  a  density  gradient  along  an  axis  normal  to 
the  output  axis,  and  if  there  is  an  applied  acceleration  component  normal 
to  both  the  gradient  axis  and  the  output  axis,  then  convection  currents  will  be 
generated  which  will  exert  a  torque  about  the  output  axis  with  resulting  drift 
rate  about  the  input  axis.  The  density  gradient  may  result  from  a  temperature 
gradient  and/or  stratification  of  a  non-homogeneous  fluid.  Most  gyro  fluids 
are  polymers  containing  a  range  of  molecular  weights  and  are  nearly,  but 
not  quite,  homogeneous. 

The  fluid  flow  must  be  determined  from  the  partial  differential 
equations  of  potential  flow  theory.  Except  for  the  very  simplest  cases,  these 
equations  must  be  solved  by  numerical  integration  methods.  There  are  many 
possible  density  gradients  occurring  in  gyros;  we  will  consider  one  simple 
example. 

Let  us  assume  that  the  gimbal  is  a  right  circular  cylinder,  and 
that  the  fluid  density  varies  linearly  along  the  spin  axis  as  expressed  b>  Eq. 
(7-1)  and  indicated  in  Fig.  7- 1(a). 

pe  *  po  [1  +  s  po  [1  +^,ln0]  (7“l) 

where 

3 

Pg  «  fluid  density  at  angle  6  -  g/cm 
P0  *  fluid  density  at  0  =  0°  -  g  /cm3 
Ap  =  P9Q  -  P2?0  -  g/cm3 
x  s  distance  along  SA  from  origin  -  cm 
R  3  radius  of  gimbal  -  cm 
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Assume  that  the  fluid  gap  h  is  very  much  smaller  than  the  radius 
R,  which  is  true  in  all  such  gyros,  and  that  the  density  of  the  fluid  is 
essentially  uniform  in  the  cross-hatched  volume  enclosed  by  the  small  angle 
AO,  gap  hf  and  gimbal  length  4 .  Let  the  gyro  be  fixed  to  the  earth  with  1A 
up,  i.  e, ,  gravity  vector  £  is  acting  downward  as  shown  in  Fig.  7- 1(a). 

Rather  than  solve  the  potential  flow  equations  which  would  require 
tedious  numerical  methods,  we  will  determine  only  the  initial  torque  exerted 
on  the  gimbal.  This  is  not  really  such  a  great  shortcoming,  since  it  is 
known  from  ocher  considerations  that  the  time  constant  for  the  convection 
flow  in  a  one-g  field  is  many  hours  or  even  several  days. 

With  the  above  conditions,  in  particular  the  narrow  fluid  gap,  the 
initial  velocity  distribution  will  be  nearly  parabolic  (approaches  parallel 
plate  conditions),  as  indicated  in  Fig.  7-l(b).  In  addition,  continuity  re¬ 
quires  that  the  initial  velocity  distribution  across  the  gap  be  the  same  all 
around  the  circumference. 

The  moment  of  the  fluid  weight  about  OA  is 
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or 


where 


2tt 

Mw  *  -Ghl  R2  ^ y  (1  +  ^-ain6)  tin  6d6 

0 


M 


-TT 


w 


■  ■  2  ■  Ghl  R  &  p 


(7-2) 


Myy  =  moment  of  fluid  weight  -  dyn*  cm 

2 

G  =  local  value  of  gravity  -  cm/s 
1  =  length  of  gimbal  -  cm 

i 

With  a  parabolic  velocity  distribution  across  the  gap,  the  vis¬ 
cous  drag  per  unit  area  on  either  the  gimbal  or  the  gyro  case  is 


4nv  6pv 

m  a 


(7-3) 


where 

2 

Fq  a  drag  per  unit  area  at  interface  -  dyn/cm 

2 

M>  =  viscosity  -  poise  a  dyn.  s/cm 

v  »  maximum  fluid  velocity  -  cm/s 
m 

v^  a  average  fluid  velocity  -  2vm/3  -  cm/s 


Formula  (7-3)  may  be  found  in  any  fluid  mechanics  textbook  for  forced  flow 
between  parallel  plates. 

The  total  viscous  drag  moment  exerted  by  the  gimbal  and  the 
gyro  case  on  the  fluid  is 

2 

a  2(2tt  RA)  FD.Ra  24\nR  A  va  (7-4) 

Because  of  the  extremely  low  fluid  velocity,  we  may  neglect 
the  inertia  of  the  fluid.  Then,  for  equilibrium  of  the  fluid,  we  have: 
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Mw  +  M*  =  0  3 


-  J  Ghl  R2  Ap  +  1  va 


or 


va  3 


Gh2  Ap 


(7-5) 


The  viscous  drag  moment  acting  on  the  gimbal  about  OA  is: 

M(lO  =  -1V2  (7-6) 


Substitute  Eqs.  (7-4)  and  (7-5)  in  (7-6). 


M  Q  =  Ghi  R2  Ap  dyn  •  cm 


(7-7) 


Notice  that  the  moment  exerted  on  the  gimbal  is  independent  of  the  viscosity 
of  the  fluid.  However,  the  time  constant  is  dependent  on  the  viscosity. 

The  resulting  drift  rate  about  the  input  axis  is: 

M  n  M  n 

u>j  =  -  -pp —  rad/sec  *  -206,  300— — deg/hr 


or 


where 


u>j  *  206.  300  Gh2  R2  Ap  =  162,  000  GhX  R2  ^ 


(7-8) 


Hf  =  rotor  angular  momentum  -  gm-  cm  /s 


Example 

Determine  the  maximum  drift  rate  of  a  gyro  fixed  to  earth  with 
OA  horizontal  and  IA  vertical  if  the  fluid  density  varies  linearly  along  the 
spin  axis.  The  gyro  has  the  following  parameters: 


=  34,000  gm»cm2/s 
h  =  0.  020  cm 
R  *1.5  cm 
i  *  2.  5  cm 


3 

Pq  *  2.  36  gm/cm 
Ap  =  0.  001  pQ 
li  *  7.  92  poises 
G  *  981  cm/s2 
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=  162,  000  x  98 1  x  0.  02  x  2.  5  x  1.  52  x 
u)j  =  i.  24  deg/hr 

Even  if  the  gyro  fluid  was  monomolecular,  which  is  not  gener¬ 
ally  true,  it  would  require  a  temperature  gradient  of  only  slightly  more  than 
one  degree  to  obtain  this  rather  large  drift  rate  in  a  one-g  field. 

The  average  velocity  of  the  fluid  in  the  gap  may  be  determined 
from  Eq.  (7-5) 


981 


x  0.  02 

-ft 


x  0.  00236  , 

x7T92 - 2,44 


10‘6  cm/ e 


With  such  an  extremely  low  velocity,  the  assumption  that  the 
fluid  inertia  may  be  neglected  is  certainly  justified  for  this  case.  It  would 
also  be  true  for  any  practical  case. 

It  is  evident  that  convection  torques  could  account  for  the  ob¬ 
served  fact  that  a  gyro  which  has  been  in  a  fixed  position  for  several  days 
(or  even  overnight)  may  take  several  hours  of  exercising  before  performance 
coefficients  are  repeatable  in  a  tumble  test.  It  would  be  interesting  to 
investigate  the  effects  of  such  density  gradients  under  the  following  conditions. 

a)  A  gyro  is  stored  at  operating  temperature  for  30  days  with 
IA  up  (as  in  a  missile)  and  then  calibrated.  Compare  the 
value  of  D<*  obtained  in  the  calibration  with  the  Dg  obtained 
when  the  gyro  is  continuously  exercised. 

b)  Repeat  the  above  test  with  SA  up  and  compare  the  value  of 
Dj  obtained  in  the  calibration  with  the  Dj  obtained  when  the 
gyro  is  continuously  exercised. 

c)  Compare  the  change  in  D~  obtained  in  (a)  with  the  change  in 
Dj  obtained  in  (b). 

If  it  is  determined  that  convection  current  effects  are  significant, 
then  it  would  be  interesting  to  determine  the  effects  on  guidance  accuracy  in 
a  high  £  environment  of  time-varying  attitude. 


APPENDIX 


GLOSSARY  OF  SYMBOLS 

In  the  general  development  of  the  gyro  equations  in  Sections  1, 

2,  and  3,  no  dimensional  units  have  been  specified  and  any  consistent  set 
may  be  used.  In  Sections  4  through  7,  commonly  used  dimensional  units 
have  been  adopted  and  it  is  those  dimensional  units  which  are  given  below. 

In  particular,  since  angular  rates  of  the  gyro  case  are  given  in  deg/hr  rather 
than  in  rad/s,  it  has  been  necessary  to  use  conversion  constants,  such  as  in 
Eqs.  (4*1)  and  (4-2),  which  are  not  required  in  the  earlier  equations. 

A  s  SI  symbol  for  ampere 

A  .  *  kth  vibratory  acceleration  amplitude  (0  to  peak)  with 

^  gyro  in  Position  No.  j  -  £ 

a  =  aT  -  £  *  applied  acceleration  vector  -  £ 

a.,  a_,  ag  3  components  of  applied  acceleration  along  I,  O,  and  S 

axes,  respectively  -  £ 

a_.,,  a0..  ,  a_  .  =  components  of  applied  acceleration  along  I,  O,  and  S 
J  ^  *  axes,  respectively,  in  gyro  Position  No.  j  and  kth 

vibration  amplitude  -  £ 

aT  3  total  acceleration  vector  with  respect  to  inertial  space  -  £ 

a  3  instantaneous  acceleration  due  to  vibration  -  g 

b  3  ratio  of  table  rate  (relative  to  earth)  to  earth  rate 

c  3  damping  coefficient 

D_  x  fixed  restraint  (acceleration-insensitive)  dx3ft  rate  • 

deg /hr 

Dt,  D^.  D-  =  linear,  acceleration- sensitive  drift  rate  coefficients  - 
1  °  5  (deg/hr)/£ 

D„,  Dcc  s  quadratic,  acceleration-sensitive  drift  rate  coefficients  - 

“  00  ss  (d.g/hr)/t2 

Dfo*  D^,  Djg  3  cross -coupled,  acceleration-sensitive  drift  rate  coef¬ 
ficients  -  (deg/hr)/£2 

EA  3  earth  axis  (polar  axis) 


Pracriiig  page  Mask 


-53- 


G 

g 

& 


& 


TT 


«r.TTr 


Hrx*  Hry*  Hrz 


viscous  drag  at  fluid -gimbal  and  fluid-gyro  case  inter* 
faces  -  dyn/cm^ 

2 

acceleration  due  to  gravity  -  cm/s 
SI  symbol  for  gram 

unit  of  acceleration  equal  in  magnitude  to  the  local  (or 
standard)  gravity  value 

gravity  vector  (in  direction  of  maximum  gravity 
gradient)  -  £ 

angular  momentum  vector  -  g.  cm^/a 

angular  momentum  vector  of  the  gimbal  with  rotor 
locked  -  g  •  cm^/s 

angular  momentum  and  angular  momentum  vector, 
respectively,  of  the  gyro  rotor  (plus  associated  rotating 
parts)  with  respect  to  the  gimbal  -  g  .  cm^/s 

components  of  17  in  the  XYZ  coordinate  system  - 

g.cm^/e  r 


h 


fluid  gap  -  cm 


I,  O,  S  s  orthogonal  coordinate  axes  along  IRA,  ORA,  and  SRA, 

respectively 

IA,  OA,  SA  2  true  input,  output,  and  spin  axes  of  gimbal,  respectively. 

Not  necessarily  mutually  orthogonal. 


IRA,  ORA,  SRA  3  orthogonal  input,  output,  and  spin  reference  axes, 

respectively,  as  indicated  by  markings  or  mounting 
surfaces  on  the  gyro  case. 


T.T.  1? 


unit  vectors  along  I,  O,  and  S  axes,  respectively 


i 


torquing  or  command-rate  current  -  A 


iONk,  iL-j.  5  command  rate-to- balance  current  in  tumble  test 

Position  No.  1,  at  table  angle  0. ,  for  counterclockwise 
and  clockwise  table  rates,  respectively  -  A 


lOSk*  iOSk 


same  as  above  except  gyro  is  in  Position  No.  2  -  A 
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command  rate -to- balance  current  in  vibration  teat 
Position  No.  j  at  the  kth  vibration  amplitude  -  A 


momenta  of  inertia  of  the  gimbal  about  the  I.  O,  and  S 
axes,  respectively  -  g.cnr 


moment  of  inertia  of  a  rigid  body  about  a  line  p  passing 
through  the  centroid 

moment  of  inertia  of  the  gyro  rotor  (plue  associated 
rotating  parts  a<Husted  for  speed  ratio)  about  its  axis  of 
rotation  -  g  •  cm^ 


principal  moments  of  inertia  of  a  rigid  body 
elastic  rsstratnt  coefficient 
command  rate  scale  factor  (deg/hr)/A 
1*  2.  3. ....  n  or  0.  1.  2. ....  n 
length  of  gimbal  -  cm 

direction  cosines  of  1.  O.  and  S  axes  with  respect  to 
the  principal  axes 

components  of  applied  normalized  torque  (applied  torque/ 
Hr)  about  the  gimbal  I  and  S  axes,  respectively 

component  of  applied  normalized  torque  about  the  gimbal 
O  axis  plus  the  normalized  viscous  drag  torque,  the 
normalized  elastic  restraint  torque,  and  the  normalized 
command  torque 

moment  o  fluid  weight  about  OA  -  dyn  •  cm 

total  viscous  drag  moment  exerted  on  the  fluid  by  the 
gimbal  and  the  gyro  case  «  dyn* cm 

viscous  drag  moment  about  OA  exerted  on  the  gimbal 
due  to  fluid  flow  -  dyn  •  cm 

integer 

radius  -  cm 

SI  symbol  for  second 


T 

3 

period  of  vibration  «  a 

TA 

3 

table  axia 

Tr  To»  Ts 

3 

component  a  of  applied  torque  about  I,  O,  and  S  axes, 
reapectively 

t 

3 

time  -  a 

T 

3 

table  rate  vector  relative  to  earth  -  deg /hr 

U,  V,  w 

3 

principal  axea  of  inertia 

U,  V,  w 

3 

orthogonal  coordinate  ayatem  fixed  to  the  gyro  caae 

u,  v,  w 

3 

vector  component  a  of  a  variable  along  the  u,  v,  and  w 
axea,  reapectively 

V 

a 

3 

average  fluid  velocity  in  gap  -  cm /a 

V 

m 

3 

maximum  fluid  velocity  in  gap  -  cm /a 

WIjk’  "Ijk 

S 

inatantaneoua  and  average  values,  reapectively,  of  the 
command  rate-to-balance  in  vibration  teat  Poaition  No. 
j  and  at  the  kth  vibration  amplitude  -  deg/hr 

WONk*  WOSk 

= 

weighted  drift  ratea  in  tumble  teat  at  table  angle  0. 
in  Poaitiona  No.  1  and  No.  2,  reapectively  -  deg/mr 

X,  Y,  Z 

- 

orthogonal  coordinate  ayatem  fixed  to  gimbal 

x.  y,  z 

3 

earth-fixed  orthogonal  coordinate  ayatem 

x.  y.  z 

3 

vector  components  of  a  variable  along  the  x,  y,  and  z 
axes,  reapectively 

o 

= 

angle  that  an  east-west  line  makes  with  a  plane  paral¬ 
lel  to  TA  and  EA  -  rad 

V  P2 

= 

angle  that  an  east-west  line  makes  with  a  plane  paral¬ 
lel  to  TA  and  ORA  at  table  angle  0,  =  0  deg  in  tumble 
test  Positions  No.  1  and  No.  2,  respectively  -  rad 

6ont»  6ost 

= 

angle  that  ORA  makes  with  TA  in  tumble  test  Positions 
No.  1  and  No.  2,  respectively  -  rad 

6te 

3 

angle  that  TA  m-h-s  with  EA  -  rad 

CI 

3 

-  <rx 
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* 

5  V 

angle  that  SA'  makes  with  the  Z  axis,  where  SA'  is 
the  projection  of  SA  on  the  YZ  plane  -  rad 

•mall  Euler  angles  required  to  rotate  the  XYZ  coor¬ 
dinate  system  into  coincidence  with  the  I  OS  coordinate 
system  when  <Pq  =  0  -  rad 

table  angle  measured  from  east-west  line  in  direction 
of  earth's  rotation  -  deg 

angle  about  OA  from  LA  to  any  point  in  fluid  gap  •  rad 
k0n,  where  k  =  1,  2,  3, , , . ,  n  -  deg 
360°/n,  where  n  is  an  integer  -  deg 

small  Euler  angles  required  to  rotate  the  UVW 
coordinate  system  into  coincidence  with  the  IOS 
coordinate  system  -  rad 

astronomic  latitude  (positive  in  northern  hemisphere)  - 
deg 

2 

viscosity  -  poise  =  dyn.s/cm 
fluid  density  at  0  =  0  deg  -  g/cm^ 

'l 

fluid  density  at  angle  0  -  g/cm 
P90  "  p270  '  S/cm3 

angular  displacement  of  the  gimbal  from  its  electrical 
null  position  -  rad 

circular  frequency  -  rad/s 

earth  rate  -  deg/hr 

angular  velocity  components  of  the  gyro  case  relative 
to  inertial  space  about  the  I,  O,  and  S  axes, 
respectively  -  deg/hr 

components  of  earth  rate  about  the  I,  O,  and  S  axes, 
respectively  -  deg/hr 

angular  velocity  of  the  rotor  relative  to  the  gimbal  - 
rad/s 
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